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Abstract 

It is well known that Pfaffian formulas for eigenvalue correlations 
are useful in the analysis of real and quaternion random matrices. 
Moreover the parametric correlations in the crossover to complex ran- 
dom matrices are evaluated in the forms of Pfaffians. In this article, 
we review the formulations and applications of Pfaffian formulas. For 
that purpose, we first present the general Pfaffian expressions in terms 
of the corresponding skew orthogonal polynomials. Then we clarify 
the relation to Eynard and Mehta's determinant formula for hermi- 
tian matrix models and explain how the evaluation is simplified in the 
cases related to the classical orthogonal polynomials. Applications of 
Pfaffian formulas to random matrix theory and other fields are also 
mentioned. 
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1 Introduction 



Pfaffian formulas were first introduced by Dyson to the theory of random 
matrices [T|. Extending Mehta and Gaudin's pioneering work [5], Dyson proved 
that the eigenvalue correlation functions for unitary random matrices were 
written in the forms of quaternion determinants, which were equivalent to 
Pfaffians. His method was then extended and applied to more general ran- 
dom matrix ensemblesj3], HJ 13 EJ |8]. Among all, Pandey and Mehta's 
work[9l Ell [EE] on the crossover between matrix symmetries revealed the 
wide applicability of Pfaffians. Extending Pandey and Mehta's result, Nagao 
and Forrester established more general Pfaffian formulas for parametric cor- 
relation functions [12]. Parametric correlation functions describe the correla- 
tions among the eigenvalues with different crossover parameters. Then Pfaf- 
fian formulas were furthermore extended to involve the multi-matrix models 
combining matrices of different symmetries |13j. In constructing Pfaffian for- 
mulas, the corresponding skew orthogonal polynomials play a crucial role. As 
a result, the skew orthogonal polynomials and related ordinary orthogonal 
polynomials are extensively used in the asymptotic analysis of the correla- 
tions. 

Recently Pfaffian formulas found new applications in combinatorics and 
related problems in non-equilibrium statistical physics. Vicious random 
walkers[m [13 US], polynuclear growth model[I7], sequences of partitions 
( Pfaffian Schur process) [18] and random involutions[19j were explored and 
statistical fluctuations were analytically evaluated. These new applications 
are regarded as discretizations of random matrix ensembles and we often 
make use of the discrete versions of orthogonal polynomials. 

In this article, we present the Pfaffian formulas employed in the analysis 
of eigenvalue correlations, list the important special cases and explain the 
relevance to random matrix theory and other applications. In §2, the gen- 
eral Pfaffian expressions are given in the forms of multiple integrals. In §3, 
we demonstrate that in a special case the Pfaffians are reduced to a simpler 
determinant expression. In §4, we explore simple cases in which the Pfaffian 
formulas become more tractable. In particular, we focus on the cases related 
to the (continuous and discrete) classical orthogonal polynomials and explic- 
itly list the corresponding skew orthogonal polynomials. In §5, the relevance 
to random matrices and other new applications are briefly summarized. 

To begin with, let us introduce the Pfaffian. For an antisymmetric N x N 
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matrix A, when N is even, the Pfaffian is defined as 

Here P = (ji, J2, • • • , Jjv) denotes a permutation of (1, 2, • • • , TV) and (— l) p 
is the sign of P. The summation J2p' is taken over all P satisfying the 
restriction j\ < j 2 , j 3 < J 4 , • • •, 3n-i < 3n- 
In terms of the vectors 

x'V' 2 *-,^ (1.2) 

with the elements 

xw = (^ ) ,4^ ) ,...,4 ) ), (i.3) 

we define the functions 

N 

p ( x« x^,...^)) = n(-r--; M) ) pfi^,^)]^,..^ 

3>l 
M-l 

x ndetb (m H^ 1 \4 m) )k«=i,2,., J v (1-4) 

m=l 

for even iV and 

AT 

P (xw,x( 2 ),...,xW) = n(4 M) -^ (M) ) 

x Pf 

M-l 

ndetb^(xf +1 \x{ m) )] i>l=1A ., JV (1.5) 



[F(xy,x\'%, l=1 , 2 ,.., N [/(4 i; )L=i,2,.,iv 

-[/ri 1) )]^l,2,.,AT 



X 



m=l 



for odd iV. We suppose that the corresponding measure is given by 

N N N 

n dM4 1} ) n <w*f ) n <w4 M) )- (i-6) 

The functions F(x,y), g( m \x,y) and f(x) are arbitrary besides the antisym- 
metry relation 

F(x,y) = -F{y,x) (1.7) 
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and the assumption that divergences are avoided in relevant calculations. As 
explained in §5, the functions p (with normalization constants multiplied) 
generalize the probability distribution functions of random matrix eigenval- 
ues. 

We are interested in the multiple integrals 

W(r {1) .-. r (1) -r {2) .-. r (2) --..-r {M) .-. r {M) ) - 1 

VV \ X 1 1 1 X k 1 > X l > -> X k 2 i , x l , , x k M ) - jjM ^ N Z h )1 

N N N 

x / n n d^(xf).../ n <w*r) 

j=fci+l j=k2+l j=kM+i 

x p(x« x( 2 \...,x( M )), (1.8) 

which give the expressions for the eigenvalue correlations. The Pfaffian for- 
mulas for these multiple integrals are presented in next section. 



2 Pfaffian Formulas for the Correlation Func- 
tions 

2.1 The case N even 

Let us recursively define 

(mn)/ / M^-?/), m = n, 

G m ' n (x,y) = j I d^^^^m-D^^)^™-!.")^^), m > n, (21) 

where S m (x) is Dirac's delta function with respect to the measure d/i m (x). 
Then we introduce 

and an antisymmetric inner product 

(f,g)^ = J d/i m (x) J dvm(y)F^ m Xx,y)f(x)g(y). (2.3) 

By means of the antisymmetric inner product, we construct the monic skew 
orthogonal polynomials R^\x) = x k + ■ • -, which satisfy the skew orthogo- 
nality relation 

(Rif\ <))W = -(<), = r 3 5 jU (2.4) 

(Rif\ Rif V M) = 0, (4fi, < } !) (M) = 0. (2.5) 



Let us 
notation 



denote the monic monomial of order j by Hj(x) = xK Then, with a 



J 3 i = {^i 



AM) 



(2.6) 



we 



u k 



find the determinant expressions for R k M \x) as 

n2fc(^) J2A; 2fc-l J2/C 2fc-2 

n2fc-l( a; ) Jjfc-l 2fc-l <^2fc-l 2fc-2 

n (x) Jo 2fe-l Jo 2fc-2 

n 2fe+ l(x) J 2 fc+1 2fe-l J2fe+1 2fc-2 ' 
n 2fc _l(x) J2fe-1 2fc-l J2fe-1 2fc-2 ' 

Hq(x) Jq 2k-l Jo 2k~2 



J2k 
^2fc-l 

Jo 



(2.7) 



i 



Jsfc+l 
Jjfc-l 

Jo 



where 



Uk 



J2k-\ 2k-l J2k-\ 2k-2 
Jlk-2 2k-\ Jjfc-2 2fc-2 



• J2/C-1 

" J2fe-2 

Jo 



+V k R { 2^\x), 

(2.8) 
(2.9) 



and Vk is an arbitrary constant 



Jo 2k-l Jo 2k-2 

. Note the indeterminacy of R^+ii 00 ) due to 



Moreover we introduce 

Rt\x) = J dn M (y)R k M \y)G^\y,x), 

S>t\x) = J d^ m (y)R ( k m) (y)F^\y,x) 
and define the matrices D^ m,n \ l( m ' n \ gO™.™) as 



(2.10) 



(2.11) 



(JV/2)-l 



1 



k=0 Tk 



(JV/2)-l 



1 



fc=0 



(m) / (™)\^f.(n.) / ("0\ (T>( Tra ) /■ ("^) \ (^1) / 



(2.12) 

-,(m,n) 



(2.13) 
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(7V/2)-l 



k=0 Tk 

where 



G 



(m,n) 



3,1 ' 

(2.14) 



F jm,n) ^M^HjW) (2.15) 

and 

r (m,n) _ / 0, m<n, 

Let us consider an antisymmetric matrix A which consists of matrix blocks 
A (m,n)^ m ^ n = i ; 2, • ■ ■ , M. Each block is a 2N x 2iV matrix which 

consists of 2 x 2 blocks 

( ) / D^™' n ^ 5*^ n,m ^ \ 

Ajj,' = c\m,n) T (m,n) ) J> ^ = 1j 2, • • • , -W. (2-17) 

Then it turns out that the Pfaffian expression for the multiple integral (11.81) 



W (JXi ... r (l). r (2) ... T (2)..... T (M) (M)s 
N/2 

= n^-l Pf ^ (m ' n) (^,^)] m ,n=l,2,..,M. (2.18) 

i=i 

Here each block y4^ m,71 )(A; m , k n ) is obtained from A*™'* 1 ) by removing the 2k m + 
1, 2k m + 2, ■ • • , 2iV-th rows and 2k n +l, 2k n + 2, • ■ ■ , 2iV-th columns. The proof 
of this Pfaffian expression is quite similar to that in [13], although it is given 
here in a slightly generalized form. 



2.2 The case N odd 

In terms of the functions G^ m,n ^ (x, y) and (x) defined in previous sub- 
section, let us employ the notations 

and 

s k = J d m (y)Ri M \y)fM(y). (2.20) 
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Moreover we introduce 



Rt\x)=Rt\x) 



-^R^\(x), k = 0,l,2,-..,N-2 

SN-l 



(2.21) 



and 



$i m) (^) 



X) 



*l m) (z)- 



Sk ^(m) 



SjV-1 



fc = 0, 1, 2, • * • , iV — 2. (2.22) 



Then the matrices Z)^ m ' n \ j( m ' n ) S^' 7 ^ are defined as 



m,n) 
hi 



j(m,n) 
1 3,l 



+ 



(AT-3)/2 



/ j \ x j l^k+lK^l ) ^k+iyXj J^k \ x i 

k=0 ' k 



(2.23) 



(7V-3)/2 



- e irW^r ] )^ Jn) 



k=0 



Tk 



X 



— K m) 1 (^ ) )/ (n) (4 n) ) - ^iri n) )/ (m) (^ m) ; 



SN-l 



(2.24) 



q(m,n) 
^3,1 



(N-3)/2 . 

= E ^ [^r ) (^ m) )^ ) + iri n) ) - ^U^)R ( 2k\x<i n) : 

k=0 ' k 



(2.25) 



As before we construct an antisymmetric matrix A so that it consists of 
matrix 
blocks 



matrix blocks A^ m,n \ m,n = 1,2, ■ ■ ■ , M and each block consists of 2 x 2 



7-)(m,n) q(n,m) 

4? n) = ( %,n) %,n) h J,/ = l,2,---,iV. (2.2G) 



It follows that the Pfaffian expression for the multiple integral (11.81) is written 

aspngun] 



W(r^ ... ^).^> ... T W..... T i JW J ... 

(JV-l)/2 

II ^-lPf[^ (m ' n) (^^n)] m ,n=l,2,.,M. 



.(1).„(2) 



„(2). 



AM) 



(2.27) 



7 



One obtains each block A^ m ' n \k m , k n ) by removing the 2k m +l, 2k m +2, • • • , 2N- 
th rows and 2k n + 1, 2k n + 2, • • • , 2iV-th columns from the block A^ m ' n \ It is 
again straightforward to prove this Pfaffian formula by following the proce- 
dures in [13] 

We remark that the above Pfaffian formulas can be extended to the case 
with the factor 

TV 

n(*r - x\ M) ) = det[(x^Y-\ l=1 , 2 ,..., N (2.28) 

(in (11 Ah and (II. 5p ) replaced by a general determinant factor [18J 

det[ Vl ^(xf ] )} u=1 ^ N . (2.29) 

In that case we need to replace rLf(a;) by <fj(x) in (12. 6p . (12.71) and (12.81) . The 
resulting skew orthogonal functions R^\x) are not in general polynomials. 



3 Reduction to Eynard-Mehta Formula 

In this section we consider the special case 

N-l i 

g^(x,y) = j-Qj{x)Rj{y), (3.1) 

where Qj{x) = x 3 ' + ■ ■ ■ and Rj(x) = x J + ■ ■ ■ are monic polynomials. The 
determinant of this function can be rewritten as 

det[</ (1) (x i ,y j )kj=i,2,-,jv 

■det[Q j _i(x/)] J - i= i i 2,...,/vdet[i? i _i(?/ z )] i j = i i 2,...,Ar 

det [(^) J-1 ]i,i=i >2 ,- 1 iv det [{yi) 3 ~\i=x,2,..,N 

N N 



1 


U N-1 


h i 


1 




U N-1 


hj 


1 





tt/V-1 l . J.J.v~J 

llj=o <h j<i j<i 



(3.2) 



so that both of the functions (11.4j) and (I1.5P are decomposed into two factors 
as 

p(x«,x( 2 )--,x( M ))=/'(xW) /V 2 \"--,x( M )). (3.3) 
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The second factor 



llj=0 % j>l 



M-1 

X 

m=2 



1 1 — ± 

J] det[^ m )(xf +1) ,x; m) )] jV=1)2 ,.., JV (3.4) 



has the form of the probability distribution function for the hermitian multi- 
matrix models. Eynard and Mehta derived a determinant formula for the 
eigenvalue correlations of such hermitian multi-matrix models[20j. Therefore, 
in this special case, their determinant formula should be reclaimed from the 
Pfaffian formulas. 

In order to see the reduction to Eynard and Mehta's formula, we choose 
Qj{x) and constants hj so that the orthogonality relation 



d m {x)d^ 2 {y)G^ M ' 2 \x,y)P j (x)Q l (y) = hjSji (3.5) 
holds. Here Pj (x) = x 3 ' + ■ ■ ■ are monic polynomials. Moreover we define 

pW( x ) = J d f i M (y)P j (y)G^ m \y,x), 1 < m < M (3.6) 

and 

QM(z) = ( d t i 2 (y)G^ 2 \x,y)Q j (y), 2 < m < M. (3.7) 



Then it can readily be seen that 

J d^ m {x) J d f , n (y)G^ n \x,y)P^\x)Q < i n \y) = h s 8 sl . (3.8) 
Let us use (13. ip to find 

G^\x,y) = fdto{z)GM(x,z)gV>(z,y) 



N-l i 

= E T-Q ( P\x)Rj{y)i 2<m<M. (3.9) 

Then we obtain 

F^ n \x lV ) = U^{x') U^{y')G^ 1 \x 1 x')G^ l \y 1 y')F{x' ) y') 



N-l N-l i 

E E 1 -rQ < f l \x)Q < i n \y)(R j ,Ri) {1 \ 2<m,n< M, 

j=0 1=0 n 3 n l 

(3.10) 
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from which it follows that 

dii M (x) J dfi M (y)F( M ' M \x,y)P^\x)P[ M \y) = (R^R^. (3.11) 

Let us choose Rk{x) so that the skew orthogonality relation 

(R 2j , R2i+i) (1) = -{Ru+i, #2i) (1) = r&u (3.12) 

(Raj, R 2 i) (l) = 0, (Raj+i, R2i + i) {1) = (3.13) 

holds. Then it can immediately be seen from (13.111) that 

Rf\x) = P ( . M \x\ (3.14) 



which yields 



R^\x) = Pj m \x), m=l,2,--.,M (3.15) 
in general. Moreover we find 

$$ m) (x) = Uii m {y)Pf>{y)F^ m \y : x) 



N-l i 

E T Q ( i m \x)(R 3 ,Ri) {1 \ 2 < m < M, (3.16) 

which means 



i=o h i 



*v\x) = T^Q%ti(x)> 2j<N-2, 2<m<M, 
mj+i 

*%Ux) = -^Q i V ) {*), 2j + l<N, 2<m<M. (3.17) 
n 2j 



Therefore, for even N, it follows from (12. 13f) and (13. 17ft that 

L J-l 



7 (m,n) = ^ 2 < m,n < M (3.18) 



and 

ift'V = -jg' m) = 0, 2<m<M. (3.19) 
Moreover, from (12.14ft and (13.171) . we obtain 

N-l 1 

S (rn,n) = £ 1 g ( m ) (x M )p (») (.,<»)) _ (jM, 2 < 7Tl,n < M (3.20) 
fc=0 ftfc 
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and 



o(m,l) 

^3,1 



0, m > 2. 



(3.21) 



We are now in a position to examine the Pfaffian expression (12.18p . It 
follows from the properties of the Pfaffian that the elements of the blocks 

A^(k m ,k n ) 



n (m,n) 
3,1 

ri(m,n) 


ri(n,m 
^1,3 




n (m,n) 
U 3,l 



<~<(n,m) 




n (m,n) 

u i,} 


°) 




' 



m, n > 2, 

m > 2, n = 1, 

m — 1, n > 2 



(3.22) 



can be replaced by 



^(m,n) 
3,^ 





c( m > n ) 



f 



' ' 



(n,m) 



, m, n > 2, 



m > 2, n = 1, 



m = 1, n > 2 



(3.23) 



without changing the value of the Pfaffian. Therefore it is straightforward to 
see that 



Pi[A^ n \k m ,k n )] m , n=1)2 ,., M 



PfL4 (1,1) (A;i,A;i)] Pf 
P^A^ik^kt)} det 





i(m,n) 



^1,3 



(n,m) 



s)y l > o 

j=l,...,fc m ,i=l,— ,fc n ,m=2,— ,M,n=2,— ,M 



Q(m,n) 
b 3,l 



j=l,— J k rn ,l=l ! —,kn,m=2,— > M,n=2,—,M 

(3.24) 



The expression ( 13.201) of S^™' 1 ^ reveals that the second determinant factor in 
the last line of the above equation gives Eynard and Mehta's formula. 
For odd N, fl3TT6|) yields 



(3.25) 
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so that 

&™\x) = $l m) (x), 2 < m < M, (3.26) 
from which we can readily find 

j{m,n) = ^ 2 < m, n < M, (3.27) 

I^V = _j(hm) = Qj 2 < m < M, (3.28) 

AT-l I 

S (rn,n) = ^ J-QM^MjpW^")) - G %> n \ 2 < m, n < M (3.29) 

fc=0 ^ 

and 

Following the steps as before, we can readily see that 



Sjft'V = 0, m > 2. (3.30) 



Pf [A^'"- 1 (k m , fc n )]m,n=l,2,-,M 

Pf^ 1 ' 1 )^!,^)] det 



j=l,-,k m ,l=l,-,k„ ! m=2,-,M,n=2,- > M ' 

(3.31) 



which gives Eynard and Mehta's formula. 

Thus we have shown that Eynard and Mehta's determinant formula in 
[20] can be derived from the special case (13 .ip of the Pfaflian formulas in §2. 



4 Simple Special Cases 

Let us go back to the general Pfaffian formulas in §2. The key ingredients of 

^( m ) 

(M), 



the Pfaflian expressions are the functions Rj^^x) generated by (I2.10p from 



the monic skew orthogonal polynomials Rj (x). We introduce the monic 

v( 

J 3 



orthogonal polynomials Cj(x) = x J ' + ■ ■ ■ with the orthogonality relation 



Now, if the function y) is expanded as 

g \X,y) 2^ (m) / ,(m+l), (m) ' ^ '> 

12 



the evaluations of the integrals in fl2.ll) and (12.101) are greatly simplified 
|2~T] . In this section, we assume this simplifying property. 

By means of the expansion (14. 2p . G^ m ' n \x,y) is evaluated as 

where we assume the completeness of C- (x) 



d m) (x)d m \y) 



Therefore, if we expand R k M \x) in terms of Cj M \ x^ as 

i=o 7j y/ij 
it can be readily seen that 

m=l,2,..;M (4.6) 

in general. Note that a>kj are independent of m. Therefore, if atkj are known, 
R ( ^ l \x) are specified for all m. It is also possible to write down the inverse 
expansion as 

u k \ x ) _ a z?M 



Y^f3 kl R\ m \x). (4.7) 

(m) Am) ] \ I V ) 

Putting (TOD into (Q leads to 

oo oo (m) (n) M™) 

Substituting H4.8[) into (12. lip and using the expansions (14. 6 j) and (|4.7p . we 
obtain 

fro 7, W 

oo oo /^( m )/,vA 

EE^^tfUi 11 )' 11 , (4.9) 
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^=0 i=u J h h 



from which it follows that 



Z=2fe+1 

<x> 



(m) 



(4.10) 



For even iV, let us put (14.101) into (12.141) . After some algebra we arrive at 



N-l (m) r {m) ( (m) , r (n) ( (n)\ 
Q (m,n) _ Tfc l x j J^fc l X l J 

- 2^ _ (n) 



AT-l oo 



(~i{rn)( ( m )\ 



fe=0 y=jV 



(m) 



(4.11) 



We also see that the expansion (14.71) can be substituted into (14. 8p . A com- 
parison with (14. 10p yields 



fc=o rfc 



(4.12) 



so that 



T(m,n) 



(m) / (m)\i(n)/ (n) 



k=N/2 



2k \ X l 



For odd N, noting 



oo oo 



(4.13) 



fc=0 i/=fc 



OO 1 



(4.14) 



we can readily find 



AT-l (m) H\n(«)/ («)\ 

ETfc u fc l X j ) U k { x l ) 
- (») 



k=0 Ik 

N-l oo 



h[ m) ht ] 



r^(m) ( ( m ) '\ 

+ E E ir v £l } ^\x\ n) ) - c^ n) 



k=0 u=N 



4^ 
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and 



M („\ (N-3)/2 



N-l 



X 



1 



E — s 2kR ( 2k + i(y) - s 2k +iR 2 n k (y) 

s n~i k=0 r k 



R { nUv) 



+ 



CO oo 



SN-1 k=N v=k 



vk 



(4.15) 



I) 



(m,n) 



E jr W^T^U^) - $£i(^ m) )^?(4 n) ) 

fc=(JV+l)/2 k 

— T~ S 2fc*2fc+lV X i J ~ s 2fc+l*2fc W J 

S iV-l fe=(AT+l)/2 rfc 

$!vn - i(y) £ 1 f^ ) (^v* + i-$a(xf ) w 



SiV-l k __ 



<N+l)/2 



r k 



(4.16) 



Thus we have shown that the Pfaffian formulas are expressed in terms 
of the expansion coefficients a 3 \ and f3jt. Let us explain how a 3 \ can be 
calculated. Suppose that R k (x) = x k + • • • are the monic skew orthogonal 
polynomials with the skew orthogonality relation 



(R 2j ,R 2l+1 ) {1) = -{R 2l+1 ,R 2j ) m =fj5 jh 
(R 2j ,R 2l )^ = 0, (R 2j+1 ,R 2l+1 )W = 

and that R k (x) are expanded in terms of Cf\x) as 



(4.17) 
(4.18) 



Rk(x) = ^a kj Cf\x), a kk = l. 

j=0 



(4.19) 



Since R k 1 \x) should equal to R k (x) multiplied by a constant c k , we can see 
that 



ttfci = akjlj^yhf^Ck, r k = c 2k c 2k+1 r k . 



(4.20) 



Noting a kk = l [ k M) \J h[ M \ we find 



Cfc = 



.,<■). 
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(4.21) 



Thus the coefficients a^j can be calculated from the expansion (I4.19p . 

Let us list the cases related to the classical orthogonal polynomials in 
which the expansions (14.191) are explicitly known. The inverse expansions 
are sometimes shown instead of (I4.19p . The inversions of the expansions are 
easy in all of the listed cases. For simplicity an abbreviated notations C n (x) 
and h n are used for C^(x) and h£\ respectively. The sign function sgn(x) 
is defined as 

{1, x > 0, 
0, x = 0, (4.22) 
-1, x<0. 

(1) Her mite [HI [S] 

When the measure is given by 



dfn(x) = I dx e~ x \ (4.23) 

J — oo 

the corresponding monic orthogonal polynomials are 

C n (x) = ±H n (x), K = ^. (4.24) 
Here H n (x) are the Hermite polynomials 

d n 

H n (x) = (-l) n e* 2 — e~* 2 . (4.25) 

The monic skew orthogonal polynomials R n (x) are given as follows. 
Case I: 

F(x,y) = e ix2+y2)/2 sgn(y - x), 

R2n{x) = C 2n (x), 
R 2 n+i{x) = C 2n +i(x) - nC 2n -i(x), 

-- 2- 2n+1 ^(2n)\. (4.26) 



Case II: 



F(x,y)=2e^ + ^/ 2 ^-5(x-y), 
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C 2n {x) = R 2n (x) - nR 2n _ 2 (x), 

C 2n +l(x) = R-2n+l, 

f n = 2- 2 "v^(2n + l)!. 



(2) Laguerre[ 

The measure is 

J dfjLi(x) — J dx x a e~ x . 
The corresponding monic orthogonal polynomials are 

C n {x) = (-l) n n\L^{x), h n = n\{n + a)\, 
where L^\x) are the Laguerre polynomials 



x a e x d n 
n! dx n 



Case I: 



Case II: 



F(x, y) = x -(*+W y -(a+i)/2 e (x+y)/2 sgn ( y _ x)? 

R 2n {x) = C 2n (x), 
R 2n +i{x) = C 2n+ i(x) - (2n)(a + 2n)C 2n -i{x), 
f n = 4(2n)!(a + 2n)!. 



F(x,y) = 2x-^/ 2 y^ a - i y 2 e^ 2 ^-5(x - y), 



C 2n (x) = R 2n {x) - (2n)(a + 2n)R 2n _ 2 (x), 
C 2n +i(x) = R 2n+ i, 

f n = (2n + l)\(a + 2n + l)\. 

An interpolation between Case I and Case II was studied in [24 
(3) Jacobi[3[6l[g[22l[23] 

The measure is 

J dm(x) = J dx (1- x) a (l + x) b . 
17 



The corresponding monic orthogonal polynomials are 

h oo+6+2n+i„| (a + n)\{b + n)\{a + b + n)\ 

n -(a + b + 2n)\(a + b + 2n + l)V 1 } 

where P^ ,b \x) are the Jacobi polynomials 

i / i \n Jn 

P^ b \x) = - — -L-J-— - {(1 - x) a+n (l + x) b+n \ . (4.35) 

n v ' (1 - x) a {l + x) b 2 n n\ dx n t K ' K > v ' 

Case I: 

F(x,y) = {l-x)-^/\l-y)-^l\l + x)-^l 2 {l + y)-^lhgTi{y-x), 



R2n{x) = C 2n (x), 
R2n+l{%) — C 2n+ i(x) 

8n(a + 2n)(b + 2n)(a + b + 2n) 
(a + b + An - l)(a + b + 4n)(a + b + An + l)(a + b + An + 2) 
_ +b+An+3( y (a + 2n)\(b + 2n)\(a + b + 2n)\ 

[ZU) - (a + b + 4n)\(a + b + 4n + 2)\ m 1 J 



C 2n -i(x), 



Case II: 

F(x,y) = 2{l-x)-^l\l-yy^\l + x)-^ 



C 2n (x) = R 2n {x) 

8n(a + 2n)(b + 2n)(a + b + 2n) 
(a + b + An + l)(a + b + An) (a + b + An - I) (a + b + An - 2) 

C 2n +l(x) = R 2n +1, 

fn = 2 a +b+ 4n +3(2n + 1} , (a + 2^ + 1)!(6 + 2n + l)!(a + b + 2n + 1) 



R 2n -2(x), 



a + b + An + l)!(a + b + An + 3)! 

(4.37) 



An interpolation between Case I and Case II was studied in 
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(4) Symmetic Hahn 

The measure is 



„ oo 

J d//i(x) = Yl 



L 



x 



L 

x 

2 



(4.38) 



where x is an integer when L is an even integer, while it is a half odd integer 
when L is odd. The corresponding monic orthogonal polynomials are 



C n (x) 



(L!) 2 (2L-2n + l)! ( - L _ w) ( L 
1 } {{L-n)\¥{2L-n+l)r n T + 2' 
n\{2L - 2n + 1)!(2L - 2n)! 



(2L-n + l)!{(L-n)!} 4 ' 
where Q^ a ' b ^(a;; L) are the Hahn polynomials |25] 

Q (a ' fe) (x- L) = J2 t^ 71 + a + b+ l )k(- x )k 1 

fc=0 



(4.39) 



{a + l) k (-L)^ 



k\ 



(4.40) 



where (a) n = (a + n — l)!/(a — 1)!. If F(x,y) is given by 

F(x,y) = 
then 

^2n(^) = C 2n {x), 

R2n+l{%) — C2n+l{ x ) /ri r . — ! \/ or ^ ; ~^2n-l( x )j 



(2L-4n + 3)(2L-4n + l) 
(2n)!(2L - 4n + 1)\{2L - An)\ 
2{2L -2n + 1)!(L -In- 1)! {(L - 2n)!} 3 ' 



(4.41) 



(5) Discrete Chebyshev 

The measure is 



r L 
/ d/ii(x) = 5^, 



(4.42) 
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where x is an integer. The corresponding monic orthogonal polynomials are 

( £ + n + !)!(„!)' 
(L-n)!(2n+l)!(2n)!' 1 ' 

where Qj£> b \x; L) are the Hahn polynomials (14.401) . If F(x, y) is given by 

F(x,y) = sgn(y - x), 

then 

C (x) - B (x) "(2n-l)(L-2n + l)(L + 2n + l) - 
C 2n (x) - R 2n (x) 2 (4n-l)(4n + l) ^n- 2 (x), 

(L + 2n + 2)!(2n)!{(2n + l)!} 3 

(L-2n-l)!(4n + 2)!(4n + 3)!' 1 ' J 

(6) Discrete exponential[19j 

The measure is 

d f i 1 (x) = Y,<l x , (4-45) 

where x is an integer. The corresponding monic orthogonal polynomials are 

q n Ti\ 

C n {x) = ,_' M n (x; l,q), 
q n (n\) 2 

hn = {l _ q) Li > (4-46) 



where M n (x; c, q) are the Meixner polynomials [25j 

M n (x;c,q) = 2 F 1 



—n, —x 

c 



1--V (4-47) 
QJ 



If F(x, y) is given by 

F(x,y) = q- {x+y)/2 a ly ~ xl/2 sgn(y - x) 
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then 



R2n(x) = C 2n (x) + g (2fc) | (l _ q)2n _ 2k C 2k (x) 



JWi(*) = C 2n+1 (x) - ^ (2n+ l)C 2n (x), 

1 - 1 - q 

(2n)!(2n + l)! v ^g 2 "+(V 2 ) 

~ (l-g) 4 " +1 (l- V^) 2 ' 

5 Relevance to Random Matrix Theory 

In this section we explain how the Pfaffian formulas are applied in random 
matrix theory and other applications. One of the most popular examples is 
the Gaussian orthogonal ensemble (GOE)[llJ. The GOE is an ensemble of 
N x N real symmetric matrices X with a Gaussian probability distribution 
function (p.d.f.) 

P(X)dX ocexp|-^TrX 2 |dX, (5.1) 

where 

N N 

dX =Y[dX j3 Y[dX 3l . (5.2) 

3=1 j<l 

Let us change the variables from the matrix elements of X to the eigen- 
values xi,X2, ■ ■ ■ ,xn and eigenvector parameters. Integrating out the eigen- 
vector parameters, we find the p.d.f. of the eigenvalues 

N N N N 

p(xi,x 2 ,- --,x N ) Y[dxj oc IJe"^) TJ \ Xj - x t \ JJ dxj. (5.3) 

3=1 3=1 3<l 3=1 



Therefore, we obtain 



N N 



p(xi,x z , ■■■,x N ) oc JJ e 2(:Ej) 11(^3 ~~ «i)Pf[sgn(ar l - a; i )]jy=i 1 2,-,jv (5.4) 

3=1 3>* 
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for even N and 



p(x 1 ,x 2 , ■ ■ -,x N ) 



N N 



[Sgn(xi - X i )]i,J=l,2,...,JV [l]i=l,2,..,7V 
- [lkl,V,iV 



OC 



iie-^n^-xopf 



i=i j>i 



(5.5) 



for odd N. These are clearly special cases of the general formulas (11.41) 
and (11.51) with M — 1. The measure and F(x,y) are given in Case I of 
(1) Hermite in previous section. Therefore the correlation functions of the 
eigenvalues 



can be evaluated in Pfaffian forms. 

The ensemble of self-dual real quaternion random matrices with a Gaus- 
sian p.d.f. is called the Gaussian symplectic ensemble(GSE). The eigenvalue 
correlations of the GSE can similarly be analyzed using the formulas given in 
Case II of (1) Hermite in previous section. It is also possible to express the 
correlation functions for a similar ensemble of antisymmetric hermitian ma- 
trices in Pfaffian forms: the derivation of the corresponding skew orthogonal 
polynomials is illustrated in [2T] . 

The eigenvalue distributions of the matrices of the form A^A, where A 
are M xN rectangular random matrices with Gaussian-distributed elements, 
are often important in applications. Here A* is the hermitian conjugate of 
A. When A are real matrices, we can employ the formulas in Case I of (2) 
Laguerre to write the correlation functions in Pfaffian forms. When A are 
real quaternion matrices, Pfaffian expression for the correlation functions are 
deduced from the formulas in Case II of (2) Laguerre. 

Hereafter we assume that iV is even for simplicity. The Pfaffian formulas 
are also useful in analyzing the dynamical correlation among the eigenvalues 
of the matrix Brownian motion model. The matrix Brownian motion model 
was also introduced by Dyson [2"o]. He considered a,n NxN hermitian random 
matrix X with the p.d.f. 




N 



Y[ dxj p{x 1 ,x 2 , ■ ■ ■ ,x N ) 
j=k+i 



(5.6) 



P(X (0) ;X;T)dX oc exp 



Tr{(X-e-X(°)) 2 } 



dX 



(5.7) 



22 



depending on the parameter r. Here X^ ^ is another matrix whose symmetry 
belongs to a subclass of the symmetry of X. Note that X is equated with 
X^> at r = 0. The measure is chosen to be 

N N 

dX = II dX n II dReXj, (UiilXj,. (5.8) 

Dyson showed that the eigenvalue p.d.f. p of X satisfies the Fokker-Planck 
equation 

where #i, #2, • • • , xjv are the eigenvalues of X and 

AT AT 

w = -Y, x ) -J2 l °g\ x j ~ x i\- ( 5 - 10 ) 

In order to solve the Fokker-Planck equation, we need to specify the initial 
condition. For example, let us suppose that X^°> is a GOE random matrix. 
Then the Fokker-Planck equation can exactly be solved and the p.d.f. p 
of the eigenvalues x^, x^ , ■ ■ ■ , xffl at r = r m , m = 1, 2, • • • , M, can be 
explicitly derived. Using the notation = (x 1 m ' ) , x 2 m \ • • • , xffl), we find 

M N 
m=l j'=l 

K ne-<4"')'/2n(,f»-xr , )Pf[F(x< i »,4 i v 1 )i^ 1A ., K 

M-l M AT 

x II det[^(^ m+1) ,x, M ; r m+1 - r m )k {=i, 2) ...,n II II ■ 

m=l m=l j=l 

(5.11) 



Here 



with 



/oo /■z' 
d*' / dz{#(x, z; -r)#(y, z'; r) - z; r')g(x, z'; r) } 
-oo J —oo 

(5.12) 



; o 
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where C n (x) and h n are defined in (14.241) and jj = j + (1/2). The dynamical 
correlation functions 



p{x ( i ) 

N 



A 



T (2). 



1 x l 1 



A 



(M)^ 



(AM) 



C=i(iV-U! 



x 



/ n y n ^f---/ n dxj%(x«,x« 



i=fci+i 



j=fc 2 + l 



(5.14) 



are then calculated in Pfaffian forms. This Pfaffian formula is known to 
be useful in analyzing the asymptotic behavior of the dynamical correlation 
functions^ [271 [28]. 

The Pfaffian formulas can also be applied to multi-matrix models in which 
one N x N real symmetric matrix is combined to N x N hermitian 
matrices X^ 2 \ X&\ • ■ • , X^^S\. The p.d.f. is given by 

P(X^,X^ 2 \---,X^)(dX) 

Af-1 



oc exp 
x exp 



-Tr lIy 1 (X«) + £ Kn(X("' 

I ^ m=2 



Tr {dlWlf 2 ' + c 2 X®X® + ■■■ + cm-iI^-^W}] (dX) 

(5.15) 

with the measure 



(dX) 



AT N 



M 



x 



AT AT 

II < II d *jf II dReXj^dlmXj 

W=2 U=l j<l 



(m) 



(5.16) 



As before we like to change the variables and integrate out the eigen- 
vector parameters. For that purpose, Harish- Chandra's integral formula can 
be utilized: for N x N diagonal matrices A and B with diagonal elements 
ai, a 2 , ■ • • , aw and b±, b 2 , ■ ■ ■ , b N , respectively, the integral over the group of 
N x N unitary matrices U can be evaluated as [29] 



J exp {— J"Tr(A - UBU ] ) 2 ^ dU 



oc 



-det 



exp 



b,) 2 



j,l=l,2,-,N 



(5.17) 
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Let us write the eigenvalues of and define 

the vectors = (xi, x 2 m \ ■ ■ ■ , x ™ )- Then, after integrating out the 
eigenvector parameters, we obtain the p.d.f. of the eigenvalues 



with 



g {m) 



P (x«,x( 2 ),---,x( a/ )) n n d 4 m) 

m=l j=l 

N 

n(xr ) -xr)pf[sgn(x / (i) -x«)],^ 1 , 2 ,.., iV 

M-l M N 

n det[^-)(xf +i \x ( r ) )\ j , l=1 , 2 ,.., N n n ^ 

m=l m=l j=l 

'x,y) = exp |--y m+1 (x) - ^V m (y) + c m xy| . (5.19) 

The correlation functions for the multi-matrix models are defined by the 
formula ( 15.141) with the p.d.f. p replaced by ( 15.181) . Comparison of the 
above formula with ( 11.41) reveals that the correlation functions are written in 
Pfaffian forms. 

Recently, through the discretization of random matrix ensembles, it turned 
out that the dynamical correlation functions for stochastic motions of inter- 
acting many particles were sometimes written in Pfaffian forms. One of 
the simplest examples is given by the vicious walkers [T4l [T5] on the lattice 
Z = {■ • ■ , —2, —1, 0, 1, 2, ■ ■ •}. Let us consider N simple symmetric random 
walkers who start from the positions 2X\ < 2X 2 < ■ ■ ■ < 2Xjy, Xj £ Z, 
at time k = 0. The positions Sj of the j-th walkers at time k > are 
restricted by the nonintersecting condition 

S{ k) <S { 2 ] <■■■ < S%\ 1 < k < K. (5.20) 

Let V(Xi, • • • , Xn; Y±, • • • , Yn) be the probability that the N walkers come 
to the positions 2Y\ < 2Y 2 < ■ ■ • < 2Y^, Yj £ Z, at an even integer time K. 
Then the nonintersecting condition implies a determinant formula 



Yj 




V(X 1 ,---,X N ;Y 1 ,---,Y N ) = —det 

(5.21) 

Hereafter we employ the scaled variables t = K/L, Xj = 2Xj/yL and r/j = 
2Yj\/L } with which a simplification occurs in the scaling limit L — > oo. 
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We consider the nonintersecting motion in the scaled time interval < 
t <T and suppose that the scaled positions of all the walkers are initially at 
the origin. Let us denote the scaled positions of the j-th walker at a scaled 
time t m by Xj and define = (x^, x 2 \ ■ ■ ■ , xffl). Then, from the 

determinant formula (15.211) . the p.d.f. of the walkers is evaluated in the limit 
L — ► oo as 

M N 

P (x<v (2) ,---,xM) n n d ^ m) 

m=l j=l 

, N N 

oc / J] dxf +1) n(4 1} - x«)Pf [sgn(x{ M+1 ) - xf +1 ^ h2 ,.., N 

J 3=1 3>l 
M M N 

n det[ g ^\x^\xi m+i % l=1 ,.., N n n d 4 m } > ^) 

m=l m=l j=l 



X 



where 



e -x 2 /(2t 1 ) e -(x-y) 2 /(2(t2-t 1 )) 



v ^rtT v /27r(t 2 - h) 

-(x-y) 2 /(2(t m+1 -t m )) 



g^ m \x,y) = 6 ] ==-, 2<m<M (5.23) 

with tjvf+i = T. The dynamical correlation functions can be defined by the 
formula ( 15. 141) with the p.d.f. p replaced by ( 15.22}) . Then we can apply the 
Pfaffian formulas in §2. The asymptotic behavior of the correlation functions 
can be analyzed by using the Pfaffian formulas |15j. 

In this article we focused on the applications of the Pfaffian to hermi- 
tian random matrices and related stochastic models. It should be remarked 
that the Pfaffian is also useful in evaluating the eigenvalue correlations of 
unitary [3D], non-self-dual real quaternion [11 j and asymmetric real[3T] ran- 
dom matrices. 
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